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Abstract 

Examples are presented of applications of the Lemaitre - Tolman model to problems 
of astrophysics and gravitation theory. They are: 1. Inferring the spatial distribution 
of matter by interpretation of observations; 2. Interaction of inhomogeneities in matter 
distribution with the CMB radiation; 3. Evolution of voids; 4. Singularities; 5. Influence of 
electromagnetic field on gravitational collapse. This review is meant to demonstrate that 
the theory is already well-prepared to meet the challenges posed by an inhomogeneous 
Universe. 

1. Why consider inhomogeneous models? 

The more appropriate question here is: why not? We know that the Friedmann - 
Lemaitre - Robertson - Walker (FLRW) models are a realistic first approximation to the 
geometry and physics of our Universe. We know (from many published examples, and see 
also below) that their generalizations are still reasonably simple. The urge to explore the 
unknown (once upon a time considered a necessary part of a scientist's personality) still 
lingers in the physics community. This should be a sufficient motivation for the beginning; 
more of it will appear along the way. 

In this article, only one class of models will be discussed, the Lemaitre (1933 and 1997) 
- Tolman (1934 and 1997) (LT) models, and only a small selection of existing material will 
be presented. For a more complete treatment see the other review by this author (Krasihski 
1997). 

2. The Lemaitre - Tolman model 

Assume that the Universe is isotropic around one observer, but not spatially homoge- 
neous. (How to save the " cosmological principle" in such a model will be explained the end 
of this section.) Assume that the matter of the Universe is dust. In comoving coordinates, 
the metric form is then: 

ds'^ = dt^ - S^{t, r)dr^ - R^it, r){d^^ + sin^ M(p'^). (0.1) 

The case R,r = has to be considered separately and it leads to the solution found by Datt 
(1938, see also Krasihski 1997). It is an inhomogeneous generalization of the Kantowski - 
Sachs (1966) solution and has interesting geometrical properties, but has not been exploited 



for astrophysics, so it will not be discussed here. Assume now R,r 7^ 0. The Einstein 
equations then imply: 

S = R,r/[l + 2E{r)]^/\ 
R,t^ = 2E{r) + 2M{r)/R + \kR^ , 

SttGp/c^ = 2M,r /{R^R,r ), (0.2) 

where E{r) and M(r) are arbitrary functions, A is the cosmological constant and p is 
the mass-density of dust. All formulae are covariant with the coordinate transformations 
r = f(r'), where / is an arbitrary function. Therefore it may be (and usually is) assumed 
that r = at the center of symmetry. The equation defining R can be solved in terms 
of the Weierstrass elliptic function (this was done by Lemaitrc already in 1933), and has 
the same algebraic form as in the Friedmann models. Indeed, the Friedmann models are 
contained in (2) as the following subcase: 

E = ^kr^, M = Mor^, R = r7^(^), (0.3) 

where k is the spatial curvature index, Mq is the Friedmann mass integral and 7?.(t) is 
the Friedmann scale factor. However, since i? is a function of both t and r, a solution of 
the second equation in (2) will depend on (t — tB{r)), where tB{r) is one more arbitrary 
function (it is a constant in the Friedmann limit). 

In the LT model every shell of constant r evolves independently of the other ones, 
with its own initial conditions, but according to the same evolution law as in a Friedmann 
model. The function M(r) specifies the initial distribution of mass, E{r) specifies the 
initial distribution of energy, and tB{r) specifies the instant of the Big Bang, which is seen 
to be position-dependent. 

The coordinate-independent definition of the Friedmann limit is: 

(^VM'),r=0 = tB,r, (0.4) 

and the functions present in eq. (4) define the two types of perturbation of the Friedmann 
model that are well-known from studies of the linearized Einstein equations: the one in 
which {E^/M'^),r 7^ = tfi^r is an increasing density perturbation, and the one in which 
{E^/M'^),r= 7^ tB,r is a decreasing density perturbation. (Increasing and decreasing 
refer to the solution in which matter expands away from the Big Bang.) 

Let us now deal with the problem of the "cosmological principle". Since the Friedmann 
models are subcases of the LT model, the functions E{r), M[r) and i_B(r) can be chosen 
so that for r > rg they go over (differentiably to a desired order) into their Friedmann 
forms. In this way, the LT model becomes an isolated island in a Friedmann background, 
and the background does not feel the presence of the island if the matching hypersurface 
r = ro is comoving. Hence, an arbitrary number of such islands in the same background 
can be considered, and such a Universe model will be spatially homogeneous in a large scale 
while being inhomogeneous in the small scale. Hence, the cosmological principle does not 
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prohibit the LT model. However, this principle is a philosophical construct and it should 
not bother anyone in exploring theoretical possibilities (see also Ellis 1984). 

This section is a very short exposition which was meant to provide a minimal back- 
ground for reading the remainig part of this article. The number of papers published on the 
LT models well exceeds 100, they are reviewed more extensively by this author elsewhere 
(Krasihski 1997). 

3. Do we really know the distribution of matter in our Universe? 

The main part of this section are results of the paper by Kurki-Suonio and Liang (1992) 
(abbreviated KSL). 

Let us assume that the LT model is used to interpret the astronomical observations. 
Observers receive their information about matter-distribution in the Universe along the 
light-cones. What they directly measure is the mass-density in the redshift space: 

^^^^ = ^''V(i)=''° ^'-'^ 

where Hq is the Hubble parameter at the observer's position at the time of observation (see 
below), z is the redshift, dV{z) = Airz'^dz is the element of volume of the redshift space 
and dm is the amount of mass in dV. (The coefficient Hq is there only for conversion to 
convenient units. The quantity p would be proportional to the physical mass-density if the 
space were flat and the Hubble law would exactly apply.) In order to infer the physical 
mass-density along the past light-cone, p{z), and the mass-density at the present time, 
p{tQ,r), two quantities have to be read out from p: the distance to the light-source, ^(2:), 
and its velocity, v{z). These cannot be measured independently, and so the problem is 
underdetermined. Without knowing I and v separately, we do not know the initial data for 
R{t, r), and so the projection from the observer's past light-cone to the hypersurface t = to 
is undetermined, too. The "knowledge" of these parameters in the Priedmann models is 
spurious - deduced from pre-assumed symmetries. 

Let us see what the LT model implies for this problem. For simplicity, it will be assumed 
that A = and that the model is a perturbation of the k < Friedmann model with the 
density parameter = 0.1. Then E > and the solution of eq. (2) can be represented 
parametrically as: 

R{t,r) = —icoshr)-l), smhri-ri = ^^^[t-tB{r)]. (0.6) 

If the observer is located at r = and the light-source at r = rgm, then the redshift and 
the Hubble parameter at the observer's position are given by: 

/■rem 

ln(l + z)= R,tr [1 + 2E{r)]-'/^dr, Ho = R,t /R 1 (0.7) 

Jo I {t,r)={to,r) 

and p is related to the physical mass-density p by: 

= ^° 4^1^'^'^'^ = 8nG il + z)z^R,r ^ ('-'^ 
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where dV = [ATrR'^R,^ /{I + 2E)^/'^]dr is the element of volume in the LT model. What 
can be said about the connection between p,p{z) and p{to,r) on the basis of information 
from observations? 

For the beginning, let us choose the increasing perturbation with the arbitrary functions 
chosen so that the function p{z) has maxima at z = 0.1 and z = 0.2 of half- width Az = 0.02 
and amplitude Ap/phom = 16 (see figures in KSL 1992), and elsewhere p/phom = 1 (^hom 
is p calculated for the Priedmann background). In this case, the physical mass-density 
along the light-cone has maxima of amplitude Ap/phom = 5.2, i.e. observation exaggerates 
reality by a factor of more than 3. Moreover, the mass-density at the hypersurface t = to 
has here maxima of unequal shape and amplitude: since the perturbation is increasing, 
and the farther condensation has had more time to evolve between the light-cone and 
the hypersurface t = to, the farther condensation has a greater amplitude and smaller 
half- width at t = to- 

As a second example, let us take a decreasing perturbation with maxima of p{z) at the 
same redshifts and the same half-width as before, but with an amplitude of Ap/phom = 
0.025. This corresponds to minima in the physical p{z) of amplitude —0.37, and to unequal 
minima in p{to,r), the farther one being broader and shallower. In this case, observation 
plainly contradicts reality. 

Since the two kinds of perturbation influence the relation between p{z) and p{z) in 
opposite ways, they can be tuned so that all inhomogencity in p{z) is cancelled in p{z). 
(The inhomogeneities in p{z) are in this case caused exclusively by "peculiar motions", not 
by inhomogeneities in the gravitational field.) Every such tuning will be unstable because 
each perturbation evolves according to a different law, and so the cancelling occurs only at 
a single light-cone; it does not survive to the hypersurface t = to, but the distribution of 
p{to,r) is different from p{z). An interesting example is given in the KSL 1992 paper: the 
function p[z) was chosen identical to the density as a function of rcdshift in the direction 
of the North Galactic Pole from the Deep Redshift Survey (Broadhurst et al. 1990), and 
the functions M{r),E{r) and tsir) could still be chosen so that p{z) was the same as in 
the unperturbed Priedmann model. 

The KSL paper 1992 concludes with a statement that the present author finds to be 
the most important memento of this whole review: ". . .it is fundamentally self-inconsistent 
to use a Friedmann relation between redshift and comoving distance when dealing with 
inhomogeneities. (...) What we see on the past light cone are only momentary "mirages"". 
This inconsistency was felt by several authors already in the 1930ies (see Krasihski 1997, 
the list would include Lemaitre, Tolman, McVittie, Dingle, W. H. McCrea and still a few 
more). The fundamentalist faith in the " cosmological principle" is of later origin. 

Considerations of similar kind were published by Moffat and Tatarski (1992 and 1995) 
and by Ribeiro (1992 and 1993). In particular, Ribeiro showed that the function p{l) (i.e. 
density vs. luminosity distance) calculated for any Priedmann model contradicts what 
astronomers expect: it is nearly constant for small distances, and decreasing for large 
distances. It is not constant because the light-cone cuts through hypersurfaces of different 
mass-densities. It is decreasing (and in fact limj^co p{l) = 0) because I — > oo when z — > oo, 
but the amount of mass within the z = oo hypersurface is finite. The conclusion from this 
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is (at A. K.'s responsibility) that if the Universe looks homogeneous to astronomers, then 
this is an observational evidence against the Priedmann models. 

4. Anisotropies in the cosmic microwave background radiation generated by 
density inhomogeneities in the Universe 

The LT model is a convenient tool for considering this problem. The inhomogeneity 
is assumed to be at and around the center of symmetry (assumed to be at r = 0). It 
may be a condensation or rarefaction, in both cases of arbitrary amplitude and density 
profile. As r tends to infinity, the functions of the LT model tend to their Priedmann 
values. Light-rays are assumed to be emitted at the time of last scattering and received 
at the present time, but they pass through the inhomogeneity along different paths (i.e. 
at different distances from the center of symmetry). The question to which an answer is 
sought here is: how does the inhomogeneity influence the temperature of the radiation 
received by the observer? The answer is found by numerical intergration of the equations 
of null geodesies in the LT model. When the tangent vector to a null geodesic is already 
calculated, the temperature at the observation event is found from: 

where "obs" refers to the observation event, "em" refers to the emission event. A;'* is the 
tangent vector to the null geodesic and is the velocity vector of the dust medium; Tem 
is the same for all rays. 

This method of investigation was first proposed by Raine and Thomas (1981) and first 
applied with a conclusive result by Panek (1992, only for the /c = Priedmann background). 
The results quoted below were obtained in the course of an elaborate project by Arnau 
et al. (1993 and 1994), Saez et al. (1994) and Pullana et al. (1996). The density ampli- 
tudes, profiles and sizes of the inhomogeneities were chosen so as to correspond to their 
current values in the observed objects (voids, galaxy clusters, the Great Attractor), but 
the r^-parametcrs of the background Priedmann regions were varied from one calculation to 
another. Any condensation in the Universe causes observer's fall toward it, and this results 
in a dipole anisotropy in the CMB radiation. This dipole and its associated quadrupole 
have been subtracted in the results quoted so that they represent the true gravitational 
perturbation. 

The anisotropy is measured by AT/(T), where (T) is the mean temperature over the 
whole sky, and AT is the difference between the actual temperature and (T) for a given 
direction. In the papers, the results are given in the form of complete graphs, from which 
only a few numbers are quoted here. Por voids, the maximal difference in AT/ (T) between 
different directions is a few times 10~^ (achieved for a model with O = 0.2, the maximal 
present density contrast —0.75, present radius of the void 30/i~^Mpc (similar to Bootes 
void) and located at redshift between 1.25 and 8.4. Por a condensation modelling the 
Virgo cluster the maximal difference is approx. 7 • 10^^ (almost identical for = 0.1 
and O = 0.2). The largest maximal difference in AT/(T), equal to approx. 4.5 • 10~^, 
was found for a model of the Great Attractor. In this last case, such models were chosen 
which produce the peculiar velocity of the observer 570 ± 60 km/s at the radial distance 
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Rq = 43^~ Mpc, while the Jl-parameter of the background, the pecuhar velocity V of 
the observer and the distance to the observer were varied. The model that yielded the 
maximal value of AT/{T) had O = 0.15, V = 570 km/s and the distance from the observer 
corresponding to the redshift 8.7. (The values of AT/{T) found by Panek for all three 
kinds of inhomogeneities, with = 1 in the background, were even lower.) The angular 
scale of the perturbation was approx. 90° the Virgo cluster, and approx. 10° for the Great 
Attractor and big voids. 

The authors of the papers quoted refrained from drawing philosophical conclusions, 
but this author is tempted to say the following. Isotropy of the CMB radiation has been 
invoked as an argument in favour of the cosmological principle already 30 years ago, when 
the precision of the measurement was at the level of 10~^. It has usually been said that 
"if the Universe were not homogeneous, then the inhomogeneities would leave an imprint 
in the CMB radiation" . Today, the error-bars stand at 10^'^' (sec Mather ct al. 1992) and 
at this level they barely make contact with the expected size of the real effect. Before 
Panek's (1992) paper, no quantitative estimates of this effect had been available and all 
such "arguments" have been uncritically repeated propaganda. Our favourite ideologies 
overshadow facts also in science, from time to time. 

5. Evolution of voids 

This application is one of the most beautiful pieces of (astro-)physics done on the 
basis of the LT model, and the so far most conclusive results have been obtained by H. 
Sato and coworkers, which justifies the presence of this contribution in this volume. The 

main references arc Occhionero et al. (1978, 1981 and 1983), Macda ct al. (1983), Sato 
et al. (1982, 1983 and 1984), Suto et al. (1984), with later contributions by Lake and 
Pim (1985), Pim and Lake (1986 and 1988), Bonnor and Chamorro (1990 and 1991) and 
Chamorro (1991). 

A void is modelled by the LT solution as follows. It is assumed that matter is distributed 
homogeneously (i.e. a Priedmann model applies) around the center of symmetry out to a 
certain distance r = ri . There it goes over into an LT region where density increases with 
distance. This region extends to r = r2, and for r > r2 again a Priedmann model applies, 
but with matter-density greater than in the central region. The boundaries between the 
three regions are not necessarily comoving with the dust. This configuration was used in the 
papers of the Occhionero and Sato groups. In the Lake approach the central homogeneous 
region contained FLRW radiation instead of dust, and in the Bonnor approach the central 
region was empty (i. e. with Minkowski geometry) and all boundaries were comoving. 

In great abbreviation, the results of H. Sato and coworkers are these: 

1. For voids, spherical shape is stable, unlike in collapse. 

2. The rim of the void propagates outward with respect to the Priedmann dust, changes 
into a high-density layer, and, in a finite time, because of pressure being zero, becomes a 
shell-crossing singularity (see next section for a definition) . Motion of the shell can then be 
studied either in the Newtonian approximation (in which case it is governed by the Sedov 
equation of blast waves) or by the Israel formalism in relativity. The singularity is not 
formed when there are nonzero gradients of pressure, but there exists no exact solution 
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which could describe this process, this case was modelled only numerically (Suto et al. 
1984). 

3. With a k = background ("background" refers to the outside Friedmann region), 
shells smaller than the particle horizon accelerate their expansion to a finite asymptotic 
velocity. Shells larger than the particle horizon let themselves be overtaken by the horizon 
and then do the same. 

4. With a A; < background, motion of the shell asymptotically freezes into the 
background flow. 

5. With a k > background, motion of the shell accelerates toward the velocity of 
light at recollapse. The void will encompass the whole space before recollapse if it was 
sufficiently large initially. 

If a void is filled with the FLRW radiation instead of a low-density dust, then it ac- 
celerates toward the particle horizon also in the A; < background (Pim and Lake 1986). 
This implies that the observed voids must be significantly younger than the Universe. 

Also here, the opposition between facts and ideology finds an illustration in the lit- 
erature. N. R. Sen (1934 and 1997) used the LT model to investigate the stability of 
the closed Friedmann model against a perturbation of matter flow which leads to a local 
decrease in density. His result, in his own words, was that "...the [Einstein and Fried- 
mann] models are unstable for initial rarefaction. . .". A similar result (instability against 
a negative-amplitude perturbation in density) was obtained by Tolman (1934). This was 
as close as one could get in 1934 to a prediction that voids should be expected to form 
in the real Universe if the Friedmann model is realistic. Unfortunately, the astronomical 
opinion-makers have known better - until voids were actually discovered in the late 1970ies. 

6. Singularities 

Two kinds of singularity can occur in the LT model: 

1. The shell-crossing singularity, at which R,r = 0. 

2. The Big Bang/Big Crunch singularity, at which R = while M,r ^ 0. 

Both are curvature singularities and p — > cxd at both. However, the shell-crossing singularity 
is considered less dangerous. It results only because p = and should go away in models 
in which pressure has spatial gradients (so far only numerical proofs of this conjecture are 
available, see Suto et al. 1984). It is not strong enough to be a problem for the cosmic 
censorship paradigm (see Joshi 1993 for a definition of strength and for a detailed discussion 
of cosmic censorship and related problems). Moreover, extensions of class through this 
singularity exist (see Clarke and O'Donnell 1993), and in the extended spacetime the shell- 
crossing singularity becomes just a caustic of dust-flow. 

In this section we shall deal with the Big Bang/Big Crunch singularity, which is un- 
avoidable when A = 0. It is immediately seen, as mentioned in sec. 2, that it is in general 
not simultaneous in the standard cosmological synchronization. The simplest consequence 
of this property is that different regions of the Universe need not be of the same age. This 
fact was exploited by Novikov (1964) and Neeman and Tauber (1967) who tried to explain 
the sources of energy in quasars as "cores of delayed expansion" - i.e. the neighbourhoods 
where the Big Bang has occurred later than elsewhere. This explanation was abandoned 
because of discrepancies between theory and observation in various details. 
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The singularities in the LT model have been investigated mostly in connection with the 
cosmic censorship conjecture. The earliest counterexample to the simplest formulation of 
this conjecture was found by Yodzis et al (1973). They showed that with appropriately 
chosen initial conditions the shell-crossing singularity can be globally naked for a finite 
period of its existence. Eardley and Smarr (1979) then found out during numerical inves- 
tigations that a globally naked singularity can be created at the center of symmetry where 
R = = M,r- With a certain set of initial conditions, in a certain moment the central 
dust particle located at i? = begins to move along a null curve on which the curvature is 
infinite. From a finite segment of this null singularity null rays can be sent to future null 
infinity, thus making the singularity globally naked. The singular null line (usually called 
"shell focusing singularity") later hits the spacelike surface of the Big Crunch singularity. 
In the comoving coordinates in which the LT model is most often represented the whole 
shell focusing singularity is mapped into a single point, making it seem to be a part of the 
Big Crunch singularity. 

Eardley and Smarr found this result for the E = LT model, later many more examples 
of such singularities were found in the E ^ models, see Christodoulou (1984), Newman 
(1986), also Krasihski (1997) and Joshi (1993). 

The Big Bang singularity in the LT model demonstrates one more kind of instability in 
the Friedmann models, first observed by Szckcrcs (1980), and then investigated by HcUaby 
and Lake (1984 and 1985b). In the Friedmann models, every light ray emitted from the 
Big Bang reaches every observer with an infinite redshift so that the Big Bang is not 
physically observable. In the LT model this is true for those locations at the Big Bang 
where tB,r ^ ^ sufficiently rapidly, and also for nonradial rays and noncomoving emitters. 
Radial rays emitted from comoving sources reach every observer with an infinite blucshift. 

The equation = is a necessary condition for a shell-crossing singularity, but not a 
sufficient one. If the set {t, rs{t)) on which R,r = has in addition the following properties: 

rs,t = = M,r{rs)=tBArs) = E,r{rs), E{rs) = 1/2, (0.10) 

then the hypersurface R,r = (which then consists of the flow lines of dust) is nonsingular, 
and it is a locus of extrema of the function R(t,r). This locus is called a neck. The 
existence of necks was first pointed out by Hellaby and Lake (1985a and c), and several 
unusual properties offered by them have been discussed by Hellaby (1987). One of these 
is the existence of a "string of beads" Universe which begins as a chain of isolated white 
holes that keep growing until they come into contact through necks, then evolve for a finite 
time, separate to form isolated black holes and collapse to their separate Big Crunches. 

As should be seen from this very short overview, also for singularities the FLRW models 
present an oversimplified picture and do not offer any insight into the more general cases. 

7. Influence of electromagnetic field on the collapse of dust 

Properties of charged dust and of neutral dust moving in an exterior electromagnetic 
field (in both cases under the assumption of spherical symmetry) have been extensively 
discussed in several papers, the most important references are Markov and Prolov (1970) 
and Vickers (1973). In the comoving coordinates the Einstein-Maxwell equations form a 
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coupled set that is difficult to handle (the set can be solved explicitly in noncomoving 
coordinates, Ori 1990). Still, it was shown that if A = and the charge density on the 
dust particles is smaller than the mass-density (in units in which c = G = 1), then the Big 
Bang/Big Crunch singularity will not occur. (Large charge density enhances gravitation 
by its contribution to the total energy-density.) The Big Bang/Crunch is avoided also for 
neutral dust with A = 0, if it moves in an exterior spherically symmetric electromagnetic 
field, created by a charge placed at the center or charges distributed in a finite volume 
around the center (Shikin 1972). In this last case the metric is given by eqs. (1) and (2), 
but the second of eqs. (2) is replaced by: 



where Q is the (constant) value of the electric charge. Hence, the electromagnetic field 
generates a repulsive gravitational field. 

This result created an expectation that a finite volume of charged dust could tunnel 
through the wormhole in the maximally extended Reissner - Nordstrom spacctimc to an- 
other asymptotically flat sheet. This expectation was killed by the papers of Ori (1990 and 
1991) who showed that while the dust is already under the inner horizon of the Reissner - 
Nordstrom spacetime, shell-crossings necessarily occur that block the passage through the 
wormhole. This happens because the dust layers that have collapsed to a smaller radius 
bounce first and collide with the next layers that still go on collapsing. The shell-crossings 
might possibly be avoided with nonzero pressure-gradients, but no such exact solutions are 
available so far. 

8. Other properties 

Several interesting properties of the LT model have been omitted in this short overview 
(see Krasihski 1997 - that one was meant to be a complete review). Among the more 

important ones are: 

1. The creation and evolution of a black hole can be studied (Barnes 1970). In the 
Schwarzschild and Kerr solutions, the black hole is stationary and present in unchanged 
form throughout the whole history of the spacetime. This does not allow for a description 
of the black hole's origin and of the process of accretion of matter onto it. In the LT 
model, one can assume such initial conditions that the black hole does not exist initially, 
then appears and keeps growing. 

2. If it is assumed that the cosmic medium (the average matter density in the Universe) 
extends throughout the solar planetary system, then circular orbits of the planets are not 
possible (Gautreau 1984). All orbits must expand by the law (calculated in the Newtonian 
approximation) : 



where H is the Hubble parameter, p is the average cosmic mass-density, M is the mass of 
the Sun and TZ is the radius of the orbit. For Saturn, this rate of expansion is 6 • 10~^^ 
meters/year - unmeasureably small, but nonzero. 




(0.11) 



^ = 87^7^^i^p/(2M), 



(0.12) 
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It is commonly believed that the last word of the relativity theory on this subject 
were the two papers by Einstein and Straus (1945 and 1946), in which a planetary system 
was assumed to be a portion of the Schwarzschild spacetime matched to a Friedmann 
background. The result was that the expansion of the Universe does not influence the 
planetary orbits. The Gautreau paper shows that this is not the only possible model. 
Moreover, the papers by Sato et al. mentioned in sec. 5 demonstrate that the configuration 
used by Einstein and Straus is unstable. The matching condition is a unique connection 
between the Schwarzschild mass and the mass removed from the Friedmann model to make 
place for the Schwarzschild region. For the k = Friedmann background the two masses 
have to be simply equal. The LT model allows one to consider the uncompensated situation, 
when the mass in the vacuole is initially larger or smaller than the Einstein- Straus value. 
This was done in some of the papers on evolution of voids - and it turned out that when the 
Einstein-Straus matching condition is not fulfilled, the edge of the void is not comoving 
with the Friedmann dust. It will either expand faster than the Friedmann background 
(when the mass in the void is "too small") or collapse toward its center. 

3. Galaxies could not have formed out of statistical fluctuations in an initially homo- 
geneous matter-density because the Universe is much too young to have allowed enough 
time for this (Bonnor 1956). In order to develop into galaxies today, the initial perturba- 
tions would have to be 10^^ times larger than the statistical fluctuations. This inference is 
usually associated with the perturbative approach to Einstein's theory which had not yet 
really existed before 1956, while Bonnor came to this conclusion on the basis of an exact 
solution. 

9. Conclusion 

The LT model was discussed here because it is the most mature one from the point 
of view of astrophysics among the several generalizations of the FLRW models that exist 
in the literature. However, many other generalizations have been found, and many papers 
were published on their properties. (The main review, Krasihski 1997, contains more than 
700 references. Readers are referred there for the literature on the subjects mentioned 
below.) In particular, the plane- and hyperbolically symmetric counterparts of the LT 
model are known, also with a nonzero cosmological constant. The same is true for the Datt 
(1938) model. All these models were generalized for electromagnetic fields of the respective 
symmetry. Generalizations with viscosity and heat-flow are also known. Szekeres (1975) 
found a generalization of the A = LT and Datt models that has no symmetry. The 
Szekeres solutions were later generalized for nonzero A and nonzero pressure (which must 
have its gradient collinear with the fluid velocity). 

The solutions mentioned above are characterized by zero acceleration and zero rotation 
of the fluid source. An even larger family of generalizations of the FLRW models exists 
in which shear and rotation are zero. One subclass of them is conformally flat, these 
are all known explicitly. In the conformally nonflat subclass, the Einstein equations were 
reduced to a single ordinary differential equation of second order; a large body of literature 
exists on particular solutions of this equation. Also in this family, generalizations with 
electromagnetic field and heat-flow are known. 

Large families of generalizations of the FLRW models are known in which 1. The 
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Vaidya radiation is present along with a perfect fluid; 2. The source is a "stiff fluid" with 
a two-dimensional commutative symmetry group; 3. The spacetime is a generahzation of 
the SenoviUa singularity-free solution. This list does not include isolated experiments with 
flnding exact solutions that cannot be organized into correlated families. 

The overall conclusion from the present article should be that investigation of inho- 
mogeneous models of the Universe is already a highly developed branch of relativity. It 
is time to begin taking all this knowledge seriously. Pretending that the FLRW models 
is everything that relativity has to say on cosmology is counterproductive and contradicts 
facts. This attitude is particularly harmful when it is spread through newly published 
textbooks. 
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